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Abstract: A campaign of 35 epochs of milli-arcsec resolution VLBA observations of the archetype colliding- 
wind WR+O star binary system WR 140 show the wind-collision region (WCR) as a bow-shaped arc of emis- 
sion that rotates as the highly eccentric orbit progresses. The observations comprise 21 epochs from the 1993- 
2001 orbit, discussed by Dougherty et al. (2005), and 14 epochs from the 2001-2009 orbit, and span orbital 
phase 0.43 to 0.95. Assuming the WCR is symmetric about the line-of-centres of the two stars and "points" at 
the WR star, this rotation shows the O star moving from SE to E of the WR star between these orbital phases. 
Using IR interferometry observations from IOTA that resolve both stellar components at phase 0.297 in con- 
junction with orbital parameters derived from radial velocity variations, the VLBA observations constrain the 
inclination of the orbit plane as 120° ± 4°, the longitude of the ascending node as 352° ± 2°, and the orbit semi- 
major axis as 9.0±0.1 mas. This leads to a distance estimate to WR 140 of 1.81 ±0.08 kpc. Further refinements 
of the orbit and distance await more IR interferometric observations of the stellar components directly. 

1 Introduction 



The 7.9-year period WR+O system WR 140 (HD 193793) is the archetype of colliding wind binary 
(CWB) systems. It is comprised of a WC7 star and an 04-5 star in a highly elliptical orbit (e ~ 0.88), 
where the stellar separation varies between ~ 2 AU at periastron to ~ 30 AU at apastron. This highly 
eccentric orbit clearly modulates the dramatic variations observed in the emission from the system, 
from X-ray energies to radio wavelengths (Williams et al. 1990). Perhaps the most dramatic variations 
are observed at radio wavelengths, where there is a slow rise from a low state close to periastron of 
a few mly, to a frequency-dependent peak in emission of 10's of mly between orbital phase 0.65 to 
0.85, before a precipitous decline immediately prior to periastron (see Fig. 1). A number of attempts 
to model these variations have been made (e.g. Williams et al. 1990; White & Becker 1995) with 
limited success, though advances in our understanding of WCRs are being made (e.g. Pittard & 
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Dougherty 2006; Pittard, these proceedings). Accurate orbital parameters are critical inputs to these 
models. 
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Figure 1: Radio emission from WR140 at 15 (blue circles), 5 (green triangles), and 1.6 GHz (red 
squares) as measured with the VLA (solid) and WSRT (open). Data taken from Williams et al. (1990) 
and White & Becker (1995). 

Many of the orbital parameters in WR 140, in particular the orbital period (P), epoch of periastron 
passage (T Q ), eccentricity (e), and the argument of periastron (u) have been established by others (see 
Marchenko et al. 2003 and references therein), and refined most recently in an extensive observing 
campaign during the 2009 periastron passage. However, the orbital inclination (i), semi-major axis 
(a), and the longitude of the ascending node (ft) require the system to be resolved into a "visual" 
binary. The two stellar components in WR 140 have been resolved using the Infra-red Optical Tele- 
scope Array (IOTA) interferometer at a single epoch (Monnier et al. 2004). This single observation 
sets the scale and orientation of the orbit since it constrains potential families of solutions for (i, a, ft). 
Further epochs of IOTA observations have been completed, but until analysis is complete, the VLBA 
observations of the WCR offer the only means to determine the orbit direction and constrain i, and 
hence ft and a. 

An initial analysis of 21 epochs of VLBA observations taken between 1999 and 2000 (orbital 
phase 0.74 to 0.95) was described in Dougherty et al. (2005). The work presented in this paper is 
an amalgamation of those observations with an additional 14 observations obtained between 2004 
to 2008, that extended the orbital phase coverage from 0.43 to 0.96. A re-analysis of the earlier 
observations with the new data leads to tighter constraints on the derived orbit parameters. 

2 VLBA observations of WR 140 

WR 140 was observed with the VLBA at 8.4 GHz at 35 epochs between orbital phase 0.43 and 0.97. 
A selection of images are shown in Fig. 2. The 8.4-GHz emission is clearly resolved, with a bow- 
shaped emission region observed at many epochs. This shape is anticipated for a WCR from model 
calculations (see, e.g., Eichler & Usov 1993; Dougherty et al. 2003; Pittard et al. 2006). The WCR 
rotates from "pointing" NNW to W over the observed orbital phases. This rotation is key to deter- 
mining the orbit of WR 140. Assuming the arc of emission is symmetric about the line-of-centres and 
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points towards the WR star, the O star is to the SSE of the WCR at epoch 0.43, and approximately to 
the E at phase 0.96. 
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Figure 2: VLBA 8.6-GHz images of WR140 at phases 0.62, 0.74, 0.78, 0.87, 0.92, and 0.93 from the 
1993-2001 orbit (phases 0.74, 0.87 and 0.93— taken from Dougherty et al. 2005) and the 2001-2009 
orbit. The synthesized beam is 2. Ox 1.5 mas 2 in the 1993-2001 orbit, and approximately 1.3x that for 
the latest observations. Note change in RA and Dec between the images taken from the 1993-2001 
orbit and the recent orbit. Rotation and proper motion of the WCR are evident during both orbits. 



3 Deriving the Orbit 

On June 17, 2003 Monnier et al. (2004) observed WR 140 to have a separation of 12.9^4 mas at a 
position angle of 151.7±1:1 degrees east of north. Using P = 2896.6 days, T Q = 2446156.3, e = 0.897 
and u = 46.8°, as determined from analysis of observations during the 2009 periastron (Fahed et al., 
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these proceedings), the observation at orbital phase 0.296 restricts potential sets of solutions for i, a 
and ft to those shown in Fig. 3 for inclinations in the range 0° < i < 180°. 
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Figure 3: Solutions for the longitude of the ascending node (red line) and orbit semi-major axis 
(green line) as a function of orbit inclination, derived from an IOTA separation and a position angle 
of the stellar components at orbit phase 0.296 (Monnier et al. 2004). The uncertainty in the IOTA 
observation gives an error in ft of closely ±1°, and in the semi-major axis of ±0.3 mas. 

The change in the orientation of the WCR with orbital phase gives the inclination since each (i, ft) 
solution family provides a unique set of position angles as a function of orbital phase. A weighted 
minimum \ 2 measure between the observed position angle of the line of symmetry of the WCR, and 
by proxy the line-of-centres of the stars, as a function of orbit phase determined for different sets of 
(i, ft) leads to a best-fit solution of % = 120° ± 4° and ft = 352° ± 2° (Fig. 4). These values lead to a 
semi-major axis of a = 8.97 ± 0.13 mas. 

4 Distance of WR 140 

Stellar distance is one of the more difficult parameters to determine for stars. Marchenko et al. (2003) 
determined a sin % = 14.10 ± 0.54 AU from radial velocity observations. From our estimate of i, 
this gives a = 16.28 ± 0.81 AU. A semi-major axis of a = 8.97 ± 0.13 mas then gives a distance 
estimate of 1.81 ± 0.09 kpc, consistent with the previous estimate of 1.85 ± 0.16 kpc by Dougherty 
et al. (2005). 

5 Summary 

High-resolution radio observations of the WCR in WR140 have provided, to date, the only way to 
determine the direction and orientation of the orbit, starting from the scale of the orbit as deduced from 
IR interferometry and orbital parameters from optical spectroscopy. From this work it is possible to 
determine a precise, and hopefully accurate, distance to WR140. 

WR140 is the best colliding-wind binary system for attempts to understand the underlying par- 
ticle acceleration processes and physics in wind-collision regions, in large part due to the wealth of 
observational constraints. However, the orbit and distance are critical to modelling the WCR, and 
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Figure 4: The change in the position angle of the axis of symmetry of the WCR as a function of orbital 
phase. The green line is the weighted best-fit curve of the position angle of the line-of-centres of the 
two stars as a function of phase, corresponding to % = 120° and ft = 352°. The other two lines show 
the fits for A\ 2 = ±1, showing the range of potential fits for the quoted uncertainties. 

further refinements of these parameters await further IR interferometric observations that will resolve 
the stellar components directly. 
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